A Short Proof of the Reducibility of Hard-Particle Cluster Integrals 
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The current article considers Mayer cluster integrals of n-dimensional hard particles in the n > 1 
dimensional flat Euclidean space. Extending results from Wertheim and Rosenfeld, we proof that 
the graphs are completely reducible into 1- and 2-point measures, with algebraic rules similar to 
Feynman diagrams in quantum field theory. The hard-particle partition function reduces then to a 
perturbatively solvable problem. 
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Introduction — Hard-particle systems are the sim- 
plest example of classical fluids developing a gaseous 
and solid phase. They are therefore a suitable starting 
point for understanding realistic fluids, as the asymp- 
totic limit to high packing densities is dominated by 
their geometry and volume dependence. Understanding 
the hard-particle free-energy structure is thus an impor- 
tant open problem. Comparing this to the better inves- 
tigated statistical and quantum field theories, one can 
identify two differences that are unique to classical parti- 
cles: their interaction potentials V(\f |) are strongly sin- 
gular, i.e. more divergent than V(r) ~ r"^ 1 in n dimen- 
sions, causing the integral J V(r)r n dr to be infinite 
over the particle's finite domain. Furthermore, the in- 
teraction is not local, resulting e.g. in the blocking of 
particles at high densities. The first problem has been 
solved by Mayer pj, by expanding the partition func- 
tion in f(rij) = exp (— PV(rij)) — 1 and representing it 
in cluster integrals. Whereas an approach to the second 
problem has been found by Rosenfeld 0-0| m splitting 
the f-function into weight functions or 1-point measures 
f( r ij) ~ M^*) ' m(^j')- His fundamental measure theory 
(FMT) of hard spheres has been extended to fluids of con- 
vex particles and to the crystalline phase [ToL [llj]. 
see [12j for a review. The weight functions entering the 
FMT, were then further investigated by Wertheim [13- 
HH, expanding ring graphs into 2-point measures. Both 
approaches make use of the observation of Kihara and 
Isihara [l7l - [i~9j that the second virial coefficient can be 
understood as the kinematic fundamental formula of inte- 
gral geometry developed by Blaschke, Santalo and Chern 
2(il - 24 1 . For further applications and an overview of inte- 
gral geometry see also [25|,[2(| and [1,0] for the discussion 
of its relation to the 1-point measures in FMT. 

In this article we will explain how the intimate rela- 
tionship between geometry and symmetry leads to the 
decoupling of arbitrary Mayer clusters into 1- and 2-point 
measures, extending Wertheim's result for ring graphs 
and justifying Rosenfeld's ansatz for the free-energy. In 
the first part, we will consider the splitting of the clus- 
ters into vertices, using their relation to integral geom- 
etry. Whereas the second part focuses onto the process 
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FIG. 1. The 4-particle cluster integral splits into two 2- and 
3-vertices along the lines A, . . . , E. 



of rejoining the vertices into 2-point measures. As the 
arguments are not restricted to 3 dimensions, we con- 
sider n > 1 dimensional particles, embedded into the flat 
Euclidean space W 1 . It will be shown that the 1-point 
measures have some similarity to the wave functions of 
bosonic field theory, where the 2-point measure can be 
interpreted as a propagator and the Mayer cluster as a 
Feynman diagram, constructed from a finite set of irre- 
ducible vertices. Each loop in the diagram introduces 
constrains that parallels charge and momentum conser- 
vation. Our discussion will focus on the mathematical 
aspects leading to the decoupling. For an extended ex- 
position with applications to 3 dimensions we refer to 

Splitting the Cluster — Mayer clusters [l| are graphi- 
cal representations of cluster integrals over the flat Eu- 
clidean space R™ with bonds as f-functions and black cir- 
cles representing space points of overlapping particles. In 
the following we will call a black circle with k outgoing 
lines a fc-vertex, corresponding to k intersecting particles 
Di fl . . . Pi Dk ^ 0. FIG. Q] gives an example for a star 
graph with 2- and 3-vertices. Each diagram is then inte- 
grated over all positions and orientations of its particles 
in R™. From an active point of view, this is identical to an 
integration over all translations and rotations in space, 
generated by the fc-fold tensor product of the n(n + l)/2 
dimensional Lie group ISO(n) ~ K™ x SO(n). Now, 
let E be a smooth, orientable Riemannian manifold of 
ri— 1 dimensions, embedded into the flat Euclidean space 
D : £ K™. Each point in the domain p £ D belongs 
to a local, orthonormal coordinate frame (p, e n , . . . , e\) 
of the tangential space. The elements of the dual space 
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then derive from Cartan's denning equations dp = CiO 1 



and de 1 



a/ e J for i,j 



1. 



The vielbeins Q % 



and spin connections w l - are the elements of the Lie al- 



gebra iso(n) with the Maurer-Cartan equations [28[ fol- 
lowing from the vanishing torsion T % — d6 l + u% A 9j — 
D(n)0i = o a nd the curvature of the flat Euclidean space 
f2y = dujij + u)^ A LOkj = 0. Here, is the covariant 
derivative in n dimensions. Let e„ be the outward nor- 
mal vector of the surface £ = dD. The curvature form 

(n) 

fiL splits then into the two contributions 



(n-l) „ 



D (n-D L 







(1) 

for a,f3 = l,...,n- 1 . (2) 



Now, the first part of the proof will show, by explicit 
calculation, that the integral density of the kinematic 
fundamental formula factorizes into an infinite set of 1- 
point measures, analogous to the known integrated form 

f i k (D 1 nD 2 n . . .) = c^ j n i (D 1 )^(D 2 n . . .) , (3) 

where the 1-point measures reduce to the finite basis fii 



of Minkowski functionals [21 1 



As a first and illustrative example of integral geometry, 
consider one particle moving in R™ . The integration over 
all rotations and translations introduces the Haar mea- 
sure of ISO(n) as the n(n + l)/2 form A™ 



<i<j 



y A "=i 



For Riemannian manifolds this differential form is iden- 
tically zero, as the vielbeins and spin connections are 
related by T % = 0. To get a non-trivial result, ob- 
serve that rotations and translations only in the nor- 
mal direction change the particle's position in space, re- 
ducing the integration to the n dimensional coset space 
SO(n)/SO(n — 1). With equation ([J), the reduced kine- 
matic measure can be written as: 



dD = A"<> Q ,„ A 9 n = 7 n ,„_ 1 iT(aD) A i 



(4) 



with the Euler class K and a constant j n ,n-i, deter- 
mined in (27j . but irrelevant for the current discussion. 
The Euler class is a topological invariant [29[ , depending 
only on the surface £ of the particle and not its domain. 
Furthermore, it is a local function, invariant under coor- 
dinate transformations of SO(n — 1). 

As already noted before, the vertex of order k corre- 
sponds to k intersecting domains D\ D . . . PI D k ^ inte- 
grated over all translations and rotations. With dD m as 
the reduced kinematic measure of particle m, the vertex 
contribution to the cluster can be rewritten 



dDi A ... A dD k | Din ... nAc#0 
= 7n,k K(d{D 1 n . . . n D k )) A dDi A ... A dD k 



(5) 



generalizing (|4]). From this we will now derive two charac- 
teristic properties of the vertices, leading to the splitting 
of diagrams into 1-point measures. 



Observe, that the Euler form of ([5]) depends only on 
the boundary of the intersecting particles. As follows 
from homology theory [23], the boundary operator is a 
derivation 9{point} = 0, d 2 = acting on homology 
elements. Its operation on two domains |22j 

d(Di n D 2 ) = dDi n D 2 + D x n dD 2 + dD x n dD 2 (6) 

can then be extended in an obvious way to an arbitrary 
set of intersecting particles. But as each operation of d 
decreases the space dimension by one, the intersection of 
n + 1 surfaces is identically zero. This proofs the first 
result that the number of irreducible terms of a vertex 
is at most n , reducing the problem to intersections of 
the form Ei n . . . n E fe n D k+ i n . . . D D m for k < n. 
This can be further simplified, as the coordinate frames 
of intersecting domains D k + 2 C\- ■ - CiD m are unrelated and 
decouple as volume forms from ([5]) . The case of £& CiD k+ i 
is similar but introduces the constrain e^ = ±ei ?£+1 ' ) 
among the normal vectors. Again, the integral measure 
of the domain factorizes from ([5]) as a volume form, with 
a sign absorbed in "f n ,k+i- It remains to determine the 
kinematic measure for £1 n ... D E^, generalizing the 
explicit calculation of Chern 22] for k = 2 to proof that 
([5]) decouples into 1-point measures. First, introduce the 
orthonormal frames (e^ a \ . . . , eff) for the a = 1, . . . , k 
particles, related by the k(k — l)/2 intersection angles 



< P ( a )| P ( b ) 



>= 



(7) 



These additional constrains generalize the coset space 
of one particle to the kn dimensional, reduced kine- 
matic measure of k intersecting surfaces ISO(n) x 
SO(fc)/ISO(n — k). In order to calculate the new n — k 
form K(Y<i CI . . . H E^), define the orthonormal frame 



(ei, . . . , e n -k,v. 



n— 5 



,v n ) 



S{e\ , . . . , e n _fc, , . . . , e,^ ^ ) , 



(8) 



where S € GL(&) is a linear, invertible transformation, 
most easily constructed by the Gram-Schmidt process. 
To parameterize the change of angular phase at the in- 
tersections, introduce the additional k(k — l)/2 angular 
coordinates < a a b < 0afc, interchanging the n — k nor- 
mal directions R{v n , v n - k +i) = (??„,..., rj n -k+i) for 
-R({o"ah}) £ SO(fe). The Euler class of k intersecting sur- 
faces is then determined by 

K{Y, X n . . . n E fc ) = A^eadth A"_ fe+1=i<j VidVj , (9) 

a n — k form in the common particle coordinates and 
the SO(k) measure of i?({av/}). But one unattractive 
feature of ([5]) still remains, as the result is not obvi- 
ously invariant under particle permutations. This can 
be solved by rewriting the connections in the principal 



,(0) 



representation u a ,h 
coordinate system of particle £1 



K^d^ and transforming to the 



After integrating out 
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a a b, one obtains a polynomial in the principal curvatures 
\K({n^}, {(f> ao })\, depending only on the intersection 
angles and the form A™=i6a ■ To get a symmetric result, 
we still have to incorporate the intersection constrains 
into the remaining integral measure dDi A ... A dDk ■ As 
explained in more detail in 27|, a suitable coordinate 
system for the domain D a follows from the above trans- 
formation (e£ , . . . , e„_ , Sen , . . . , Sen ), correspond- 
ing to the decomposition SO(fc) = SO(fc)/SO(fc — 1) x 
. . . xSO(3)/SO(2) x SO(2). Taking into account the ad- 
ditional constrains 6^ = at E a , the integral measure 
for the k- vertex obtains the final form 



dDi A ... A dDk | El 



n...nS fc7 to 



= J n ,k\K{{K^}, {</> a6 })| dG A dSi A ... A dS fe 



(10) 



factorizing into the reduced kinematic measures of in- 
dividual particles dH m , the group measure dG and the 
determinant \K\. 

Equation (|10| generalizes Wertheim's decoupling of the 
second virial diagram [13] into 1-point measures p, l (D a ) 
by expanding the (j} ao dependent function |if|dG into ten- 
sor products of its vectors. Formally, this splitting par- 
allels the known generalization of the integrated form of 
the kinematic fundamental equation of integral geometry 
(131) . But the Minkowski functional [2l[ are now replaced 
by the infinite set of 1-point measures 



4(D)=p (T ({K a })l[(e i ) 



® 7T 



for the f- function A at a fixed root point (see FIG. [TJ, 
where the multi-index k — (a, M) carries further infor- 
mation about the polynomial p a of principal curvatures 
K a , classified by the Young diagram a of the symmetric 
group and the tensorial exponents Mj £ N* of the vec- 
tors. The additional variable n indicates the parity of 
the tensor valued function under axial rotations, defined 
later. 

This infinite tensor space is the result of an explicit 
manipulation of (|10p . But there is one alternative ap- 
proach that might hint at connections to index theory 
29]. From the Gauss-Codazzi equation (|2|) we see, that 
K a , e a and e n are covariantly constant and lie in the ker- 
nel of the eigenvalue equation _D(™ _1 )$ = X<&. Therefore, 
we can characterize the set of 1-point measures as the 
tensor valued space 

M(E, M") = {$(k q , <x) I D^-VQ = 0} (12) 

leading to the second result that the irreducible elements 
of k-vertices decompose into the infinite dimensional ten- 
sor space of 1-point measures yV((E,R") with multiplica- 
tive structure defined by the kinematic equation H0\) . 
Comparing this to the Hilbert space of quantum mechan- 
ics, one might interpret the 1-point measures as the wave 



functions 4> of a massless bosonic particle with vertex in- 
teractions up to order <f> n . 

There are many open questions concerning this space, 
of which some will be discussed in 27[ . Let us here only 
mention the case of Riemannian surfaces T g of genus g 
and Euler characteristic x{T g ) = 2 — 2g. For g > 1 the 
volume dependent part of the second virial x(-Di) v ol(-D2) 
is either zero or negative, contradicting experience. This 
inconsistency follows from the representation in Mayer 
functions and can be mended by including further ho- 
motopic invariants as the Euler linking number [29j |. As 
these forms are again decomposable into 1-point mea- 
sures, the above result remains unchanged even for ho- 
motopically nontrivial manifolds. 

Rebuilding the Cluster — With the decoupling of the 
cluster integrals into vertices, diagrams as that of FIG. Q] 
reduce to expressions of the form 



r n[AE r ED]h r 
^iih Sifci °J2»2 



AB^nlBE r EC]l 3 r n CD 
S'sJa ^ksis ^jiU 
A ,.A\I ,.B ,.B\( ,.C ,,C \i ,,D ,.D\i ,.E ,.E 



(13) 



where the angular brackets between the coefficients C,^ 
indicate the symmetric permutation of f-labels. Calcu- 
lating the cluster integral is now reduced to integrat- 
ing pair products of 1-point measures {p^pf^ for the 
rigid particle S^. But both measures have been de- 
rived at independent points in the embedding space, 
here denoted by f\ and r-i- It is therefore necessary 
to merge these two coordinate systems into the one of 

e„), with the axial 



(11) the rigid body, denoted by (ei, 



direction z n r = f\ — r%. Generalizing Wertheim's dis- 
cussion to n dimensions, the tensorial part of (llip can 
be reduced by integral averaging over axial rotations, re- 
sulting in breaking up the symmetry group SO(n) into 
SO(n)/SO(n - 1) x SO(n - 1) = S"" 1 x SO(n - 1). But 
the two Lie groups SO(2Z + l) = B l and SO(2Z) = D t be- 
long to different Dynkin diagrams, where Di carries an 
additional Z2 automorphism [28| that induces possible 
sign changes under axial rotations. This is an additional 
degree of freedom, we have taken care of by including 
the parity symbol n in (|11|) . The existence of such odd 
and even parity tensors might have a significant influence 
on the solid phase structure, as the free-energy functional 
could develop further minima in addition to the pointwise 
vanishing of the group measures in (fT0|) by sin (4> a b) = 0. 
It would therefore be interesting to compare the odd to 
the even dimensional phase structure, with a special fo- 
cus on n = 8 and its Z3 automorphism group of D4. 

Calculating the axial average of products of SO(n) in- 
variant tensors A, mi, . . . , mi >, reduces now to deter- 
mining their SO(n — 1) invariant subspaces 



< X,rh'\B\X,m > = const, for BeSO(n-l) 
A(m',m) = 



(14) 



represented by linear constrains between their weight vec- 
tors A(m', to) = 0. A general solution of this problem has 



4 



been given by Cartan [28| and extended to a complete 
classification of Riemannian symmetric spaces. We can 
therefore assume the constrains A = to be known and 
indicate axially symmetric tensors of rank A and weight 



vector m by |A, rh, e > 



A=0> 



where the third component 



indicates either the coordinate system of the vertex e, the 
rigid body e or that of the embedding space E. Following 
the notation of 13Hl6j. the 2-point measures are tensor 
products of the form 

Kn) ® v(fj) = 

f{r lJ ,X l ,X 3 ,M l ,M J )\X i ,M i ,e l > ®\\j,Mj,ej > , 

(15) 

whose tensorial part will also be written as |e, > ®\ej > 
if no further reference to the representation space is re- 
quired. Joining the two coordinate systems ei,ej into 
that of the rigid body e, corresponds to an integration 



J < t\ei >< t\ej > d>jj l dbjj |e > <8)|e > 



(16) 



e > ® e > U=o 



over the Haar measure duj of SO(?i — 1), reducing the 
space by (n — l)(n — 2)/2 dimensions. The resulting 2- 
point measures 



/(r)|e > <g>|e > | A =o 



(17) 



are now irreducible tensors in the coordinates (e, r) 6 



S" 



of the particle S. 



The 2-point measures and the expansion coefficient of 
the vertex ([3]) carry the local properties of the particle S. 
But the cluster diagrams also carry global informations, 
as can be seen from FIG. [TJ Consider the star graphs of 
the free-energy. With the exception of the second virial, 
all such diagrams are build from closed subdiagrams or 
loops, introducing a set of constrains r\2 + . . . + r m i = 
between the coordinates of the m root points of one loop. 
These constrains also couple particles that otherwise do 
not intersect. This problem occurs first for ring graphs 
and has been solved by Wertheim in 13l4l6j by introduc- 
ing the Radon transformation [30| to replace the parti- 
cle fixed coordinate system (e,r) by that of the embed- 
ding space (E,R). This approach can be extended in a 
straightforward way to arbitrary diagrams in n dimen- 
sions. The Radon transformed 2-point function 



F(R,E) 



J < E\t >< E\t > \ A= of(r)5(rcE - R)d n x 



(18) 



can then be seen as the propagator of a state vector \E > 
along a line of distance R. With \E > F(R, E) < E\ in- 
terpreted as the Green's function, the Mayer cluster can 
be seen as the analog of Feynman diagrams, with the 



constrains 

(r 12 ei2 + ...+r m , 1 e„ l ,i)£; = 0, A = (19) 
as the conservation laws for each loop. Calculating the 
virial coefficient reduces now to an integration over the 
embedding space coordinates (R,E), completing the ar- 
gument that arbitrary hard-particle clusters in n > 1 
dimensions decouple into 2-point measures. The advan- 
tage of this approach is not an efficient calculational tool 
but the detailed connection between the 1-point mea- 
sures and their occurrence in the free-energy, opening up 
a path to better understanding the phase structure of 
hard-particle fluids by the particle geometry and relative 
orientations. 
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